Search for extraterrestrial point sources of neutrinos with AMANDA-II 
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We present the results of a search for point sources of high energy neutrinos in the northern 
hemisphere using AMANDA-II data collected in the year 2000. Included are flux limits on several 
AGN blazars, microquasars, magnetars and other candidate neutrino sources. A search for excesses 
above a random background of cosmic-ray-induced atmospheric neutrinos and misreconstructed 
downgoing cosmic-ray muons reveals no statistically significant neutrino point sources. We show 
that AMANDA-II has achieved the sensitivity required to probe known TeV 7-ray sources such as 
the blazar Markarian 501 in its 1997 flaring state at a level where neutrino and 7-ray fluxes are 
equal. 



INTRODUCTION 

The search for sources of high-energy extraterrestrial 
neutrinos is the primary mission of the Antarctic Muon 
and Neutrino Detector Array (AMANDA). The mecha- 
nism for accelerating cosmic rays to energies above the 
"knee" (10 15 eV) remains a mystery. Cosmic rays are 
thought to be accelerated in the shock fronts of galactic 
objects like supernova remnants, microquasars and mag- 



netars, and in extragalactic sources such as the cores of 
active galaxies (AGN) and gamma ray bursts (GRB) 0. 

High energy protons accelerated in these objects will 
collide with the ambient gas and radiation surrounding 
the acceleration region, or with matter or radiation inter- 
vening between the source and the Earth. This leads to 
pion production, the charged pions decaying into highly 
energetic muon and electron neutrinos, and the neutral 
pions decaying into the observed 7-rays. Fermi accel- 
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eration of charged particles in magnetic shocks naturally 
leads to power-law spectra, E a , where a is typically close 
to -2. By the time the neutrinos reach the Earth, vacuum 
oscillations will have uniformly populated all three flavors 
(unless the neutrinos are unstable f]). All limits quoted 
in this letter are on the muon neutrino flux arriving at 
the Earth; limits at the source will be approximately a 
factor of two higher due to oscillations. 

Neutrino astronomy may provide information comple- 
mentary to the knowledge obtained from high energy 
photons and charged particles, since neutrinos propagate 
directly from their point of origin undeflected by mag- 
netic fields. Therefore they have the potential to reveal 
"hidden" sources masked by photon absorption. Probing 
the neutrino sky may bring us closer to solving the cosmic 
ray mystery, or might even reveal something completely 
new and unexpected. 



AMANDA-II 

AMANDA-II ;J| is a Cherenkov detector frozen into 
the antarctic polar icecap. A high energy muon neutrino 
interacting with the ice or bedrock in the vicinity of the 
detector results in a high energy muon propagating up 
to tens of kilometers. 

The muon track is reconstructed based on detection of 
the Cherenkov light emitted as it propagates through a 
19-string array of 677 photomultiplier tubes. The median 
neutrino pointing resolution is 2 - 2.5°, depending weakly 
on declination, and is dominated by the resolution of the 
muon track reconstruction. AMANDA-II demonstrates 
a significant improvement over its predecessor in accep- 
tance and background rejection, especially near the hori- 
zon. Results from the first phase of AMANDA, the 10- 
string subdetector AMANDA-BIO, have been reported 
in 0. 

Atmospheric muons from cosmic rays that penetrate 
to AMANDA depths are the dominant background. 
AMANDA-II views the neutrino sky above the north- 
ern hemisphere using the Earth as an atmospheric muon 
filter. Cosmic rays also produce neutrinos in the Earth's 
atmosphere, but with a spectral index of a = —3.7, 
softer than expected for astrophysical sources. Atmo- 
spheric neutrinos are an important source for calibration 
in AMANDA , but are also background to a search for 
extraterrestrial point sources. A point source search is 
conducted by looking for excess events above the back- 
ground, which is experimentally measured for a given an- 
gular search bin by taking the average background rate 
in the same band of declination. 



DATA PROCESSING & DETECTOR 
SIMULATION 

The data set comprises 1.2 x 10 9 triggered events col- 
lected over 238 days between February and November, 
2000, with 197.0 days livetime after correcting for 17.2% 
detector dead-time. After application of an iterative 
series of maximum-likelihood reconstruction algorithms, 
2.1 million events reconstructed with declination 6 > 0° 
remain in the experimental sample. 

To prevent bias in the selection of cuts, the data are 
"blinded" by randomizing the reconstructed Right As- 
cension (R.A.) angle of each event. The output of a Neu- 
ral Network (NN) trained on simulated events and using 
six input variables (such as the number of unscattered 
photon hits, track length, likelihood of the muon track 
reconstruction, and topological variables 0), is used as 
a quality cut. A second cut is placed on the likelihood 
ratio (LR) between the muon track reconstruction and a 
muon reconstruction weighted by an atmospheric muon 
prior 0- The prior describes the zenith-dependent fre- 
quency of downgoing muons such that choosing a cut 
on the LR gives downgoing hypotheses a prior weight of 
up to 10 6 more than the upgoing hypotheses, effectively 
forcing events surviving the cut to be of higher quality. 
The final choice of NN quality cut, likelihood ratio cut, 
and the optimum size for a search bin are determined 
independently in each 5° band of declination in order 
to optimize the limit setting potential of the experiment 
|5j. The directional information is then restored (data 
"unblinded" ) for the calculation of the limits and signif- 
icances. 

A full simulation chain 3] including neutrino absorp- 
tion in the Earth, neutral current regeneration, muon 
propagation and detector response is used to simulate 
the point source signal according to an E~ 2 energy spec- 
trum. The limits obtained in this analysis are a func- 
tion of the measured background, rib, as well as the ex- 
pected number of events, n s , from a simulated flux Q(E): 
^linut(-E) = $(E) x [i 90 (n obs ,n b )/n s where n bs is the 
number of observed events in the given source bin, and 
^90 is the 90% upper limit on the number of events fol- 
lowing the unified ordering prescription of Feldman and 
Cousins @. 

SYSTEMATIC UNCERTAINTIES 

Atmospheric neutrinos were used to determine the ab- 
solute normalization of the detector simulation. The 
normalization factor 0.86 is consistent with the theoreti- 
cal uncertainty of 25% on the atmospheric neutrino flux 
4J. The overall systematic uncertainty, which includes 
the theoretical uncertainty of the atmospheric neutrino 
flux and statistical uncertainty of the measured back- 
ground, is incorporated into the limits using the Cousins- 
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Highland 191 prescription with unified Feldman-Cousins 
ordering |8lll0| but with a more appropriate choice of the 
likelihood ratio test 

Coincident events between the SPASE air shower array 
[l2| and AMANDA-II were used to evaluate the system- 
atic error in pointing accuracy. This value was deter- 
mined to be less than 1°, which results in signal loss in a 
typical search bin of only 5%. 



RESULTS 

The final sample consists of 699 upwardly recon- 
structed (8 > 0°) events, illustrated in Fig. ^ A compar- 
ison to the normalized atmospheric neutrino simulation 
reveals that for declinations S > 5° the sample is strongly 
dominated by atmospheric neutrinos. 
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FIG. 1: Final point source search sample plotted in equatorial 
coordinates. The thick band of events at S < shows the 
onset of cosmic muon background contamination. 

A binned search for excesses in the region 0° < S < 85° 
has been performed. The search grid contains 301 rect- 
angular bins with zenith-dependent widths ranging from 
6° to 10°, based on the aforementioned bin-size optimiza- 
tion. The grid is shifted 4 times in declination and 
right ascension to fully cover boundaries between the bins 
of the original configuration. The number of times to 
shift the grid was studied, taking into account statisti- 
cal penalties for the number of shifts by using simulated 
event samples, and set at a level where further shifts do 
not markedly improve the average maximum significance 
obtained on simulated Poisson-fluctuated signals of sim- 
ilar magnitude to the background. 

The most significant excess, observed at about 
68° Dec, 21.1hR.A., is 8 events observed on an expected 
background of 2.1. Simulation reveals a probability of 
51% to observe such an excess as a random upward fluc- 
tuation of the background. 



In addition to the binned search, we place limits on a 
number of extragalactic and galactic candidate sources. 
Circular bins with optimized radii are positioned at each 
candidate position; the number of expected background 
events is given by the number of observed events in the 
declination band scaled down to the search bin area. The 
same method applied to any point in the Northern hemi- 
sphere yields neutrino flux upper limits shown in Fig. 
For the region S > 85° an adjacent region at lower dec- 
lination was used for the background estimation. The 
average flux upper limits for E~ 2 spectra obtained in an 
ensemble of identical experiments in the case of no true 
signal is shown vs. declination in Fig. It should be 
pointed out that due to the large number of cells tested 
a single source with a flux at the sensitivity level (90% 
average flux upper limit) would normally be interpreted 
as a statistical fluctuation. Therefore a flux a few times 
higher would be needed to make a discovery of a point 
source possible in the binned search. 




a [h] 



FIG. 2: Neutrino flux upper limits (90% confidence level) in 
equatorial coordinates. Limits (scale on right axis) are in 
units of 10 _7 cm~ 2 s _1 for an assumed E~ 2 spectrum, inte- 
grated above E v — 10 GeV. Systematic uncertainties are not 
included. 

In Tab. [IJ we present neutrino flux limits for north- 
ern hemisphere TeV blazars, selected GeV blazars, mi- 
croquasars, magnetars, and selected miscellaneous can- 
didates. The limits are computed based on an assumed 
E~ 2 energy spectrum. Limits for other spectra can be 
computed using the neutrino effective area shown in Fig. 
0] The neutrino flux limit for an assumed flux d<I> / dE oc 
E a is inversely proportional to the energy averaged effec- 
tive area Atf{a) = J^A cS {E)E a dE/ J^E a dE, i.e. 
a limit for spectral index a ^ — 2 is obtained by multi- 
plying the presented limit by A c s (~~2)/A c g (a). Effective 
areas at declinations not shown in Fig. Q]can be obtained 
by linear interpolation in 5; the systematic shift induced 
by this interpolation is below 20% for spectral indices in 
the range a = — 1.5. . . —2.5. 
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FIG. 3: Average flux upper limits (90% and 99% confidence) 
for an E~ 2 signal hypothesis, integrated above E v — 10 GeV, 
are shown vs. declination (<5), solid lines with, and dashed 
lines without the inclusion of systematic uncertainty (as de- 
scribed in the text). The limits worsen near the horizon due 
to the onset of cosmic ray muon contamination. 



In [13(|, expected numbers of neutrino induced muon 
events for galactic microquasars using source parame- 
ters estimated from radio observations are calculated. 
For microquasars whose jets have not been resolved in 
the radio band the neutrino emission is estimated from 
the synchrotron luminosity. In the case of the micro- 
quasar SS433, 252 muons yr _1 km~ 2 are predicted. Scal- 
ing to the AMANDA-II effective area at that declination 
(Ag ff = 7900 m 2 ) and to the livetime of this analysis, 
yields a prediction of 1.07 events for an assumed E~ 2 
spectrum. We observe no events in the search bin for 
this source, and place a 90% upper limit of 1.24 events. 
Due to a random fluctuation, this limit is about 3 times 
better than the sensitivity at this declination. 

In Fig. [S] the AMANDA-II neutrino sensitivity is com- 
pared to the 1997-averaged TeV 7-ray flux of the blazar 
Markarian 501 (z — 0.031), and the intrinsic source spec- 
trum (corrected for IR absorption). The figure demon- 
strates AMANDA-II has achieved the sensitivity needed 
to search for neutrino fluxes from TeV 7-ray sources of 
similar strength to the intrinsic 7-ray flux. 

Data collected in 2001-2002 are being analyzed and 
should improve the sensitivity of this analysis approxi- 
mately by a factor 2.3. 



( QTinln dtfl 
V_y cl 1 1 LI 1 LI a L e 


Dec. [°] 


R.A. [h] 


^obs 


lib 


<3E> lim 


(T> lim 




TeV Blazara 










A/TarVnrian 491 


38.2 


11.07 


Q 
O 


1 50 


3 


'\ 5 


iYTfiT*i<"ftT"i£iTi i~ini 

1.VJ. CLl XYCI.L ±CL11 UU1 


39.8 


16.90 


1 


1.57 


1.5 


1.8 


1 ES 1426+428 


42.7 


14.48 




1.62 


1.4 


1.7 


1 P,S 2344+51 4 


51.7 


23.78 


1 


1.23 


1.6 


2.0 


1 ER 1 959+650 


65.1 


20.00 





0.93 


0.9 


1.3 




Ge V Blazars 










OSO 0^98-1-1 *34 


13.4 


5.52 


I 


1 HQ 


9 5 

Z . O 


9 n 

Z--U 




16.6 


2.62 


I 


1.49 


9 n 


1 7 


OSO 1 fi1 1 -1-^4^ 


34.4 


16.24 


n 
u 


1.29 


n 7 


n 8 

u.o 


OSO 1633+382 


38.2 


16.59 


1 


1.50 


1.5 


1.7 


ORO 0219+428 


42.9 


2.38 


x 


1.63 


1.4 


1.6 


ORO 0954+556 


55.0 


9.87 


I 


1.66 


1.3 


1.7 


ORO 0716+714 


71.3 


7.36 


2 


0.74 


2.9 


4.4 




Microquasars 












5.0 


19.20 


n 
u 


2 "38 


1 n 


n 7 




10.9 


19.25 




n qi 


9 Q 
z. ,7 


2.2 


HRD 70499-1-^9 


32.9 


4.36 


2 


1 "^1 


9 Q 
z. .7 


9 

Z . 3 




35.2 


19.97 


9 

Z 


1 Id 


9 9 
z. z 


9 5 
z . J 


Cygirus X3 


41.0 


20.54 


3 


1.69 


3.0 


3.5 




48.0 


11.30 


1 


0.92 


1.7 


2.2 


OT Cam 

V. . 1 V 'dill 


56.0 


4.33 





1.72 


0.6 


0.8 


LR I +61 303 


61.2 


2.68 





0.75 


1.0 


1.5 


01 ViX 


magnetars 


& miscell 


aneous 






rpt? ionn+14 


9.3 


19.12 





0.97 


1 4 

1.4; 


i n 

1 .u 


f~V»tt NToKnla 


22.0 


5.58 


2 


1.76 


9 6 
z.u 


9 4 

Z .4: 


vycihDiupeia ^1 


58.8 


23.39 





1.01 


n q 


1.2 


SF,0 Tfl45n+11fl5 


11.4 


4.82 


2 


0.89 


4.2 


'1 9 
■ z 


M 87 


12.4 


12.51 





0.95 


1 '-t 


i n 


Geminga 


17.9 


6.57 


3 


1.78 


3.7 


3.3 


UHE CR Triplet 


20.4 


1.28 


2 


1.84 


2.4 


2.3 


NGC 1275 


41.5 


3.33 


1 


1.72 


1.4 


1.6 


Cyg. OB2 region. 


[15J 41.5 


20.54 


3 


1.72 


2.9 


3.5 


UHE CR Triplet 


56.9 


12.32 


1 


1.48 


1.4 


1.9 



TABLE I: 90% upper limits on candidate sources. The num- 
ber of events observed within the search bin is denoted by 
n Q ba, and rif, is the number of expected background events 
determined by measuring the background off-source in the 
same declination band. Limits are for an assumed E„ 2 spec- 
tral shape, integrated above E v — 10 GeV, and presented in 
units of 10~ 15 cm~ 2 s^ 1 ($ M ) and lO^cm - ^" 1 ($„). 
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FIG. 4: Neutrino and muon effective areas vs energy at different declinations (S). is the muon energy at closest approach 
to the center of the detector. The effect of neutrino absorption in the Earth is included in the neutrino effective areas. 
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FIG. 5: The AMANDA-II 90% average flux upper limit (197 
days livetime) for two assumed spectral indices (a) is com- 
pared to the average 7-ray flux of Markarian 501 as observed 
in 1997 by the HEGRA system of air Cherenkov telescopes 
fl4|. These average upper limits are based on the assump- 
tion that the neutrino spectrum extends to beyond 10 PeV. 
Also shown is the intrinsic source flux after correction for IR 
absorption by de Jager and Stecker [Ttijl . The shaded area 
is bounded by two curves corresponding to different mod- 
els of galactic luminosity evolution. For comparison, the 
AMANDA-B10 result is also shown. 
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